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ABSTRACT 

Aims. To explore the nature of the unidentified very-high-energy (VHE, E>100 GeV) gamma-ray source HESS J1626-490, we 
investigated the region in X-ray, sub-millimeter, and infrared energy bands. 

Methods. So far only detected with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes, HESS J1626-490 could not 
be unambiguously identified with any source seen at lower energies. Therefore, we analyzed data from an archival XMM-Newton 
observation, pointed towards HESS J 1626-490, to classify detected X-ray point sources according to their spectral properties and 
their near-infrared counterparts from the 2MASS catalog. Furthermore, we characterized in detail the diffuse X-ray emission from a 
region compatible with the extended VHE signal. To characterize the interstellar medium surrounding HESS Jl 626-490 we analyzed 
12 CO(J=1-0) molecular line data from the Nanten Galactic plane survey, HI data from the Southern Galactic Plane Survey (SGPS) 
and Spitzer data from the GLIMPSE and MIPS GAL surveys. 

Results. None of the detected X-ray point sources fulfills the energy requirements to be considered as the synchrotron radiation 
counterpart to the VHE source assuming an inverse-Compton (IC) emission scenario. We did not detect any diffuse X-ray excess 
emission originating in the region around HESS Jl 626-490 above the Galactic background and the derived upper limit for the total 
X-ray flux disfavors a purely leptonic emission scenario for HESS Jl 626-490. We found a good morphological match between 
molecular and atomic gas in the -27km/s to -18km/s line-of-sight velocity range and HESS J1626-490. The cloud has a mass of 
1.8xl0 4 M o and is located at a mean kinematic distance of d = 1.8 kpc. Furthermore, we found a density depression in the HI gas 
at a similar distance, which is spatially consistent with the SNRG335. 2+00.1. We discuss various scenarios for the VHE emission, 
including the CO molecular cloud being a passive target for cosmic ray protons accelerated by the nearby SNRG335.2+00.1. 

Key words. Acceleration of particles - ISM: supernova remnants - ISM: clouds - ISM: individual objects: HESS Jl 626-490 - X-rays: 
ISM - Submillimeter: ISM 



1. Introduction 

During scans of the Galactic plane with the Imaging 
Atmospheric Cherenkov Telescopes (IACTs) of the H.E.S.S. 
(High Energy Stereoscopic System) array a number of new 
Galactic very high energy (VHE, E>100 GeV ) y-ray sources 
were detected (lAharonian et al.ll2005ai l2006al) . Many of these 
sources could be identified as pulsar wind nebulae (PWN), 
shell-type supernova remnants (SNR), y-ray binaries, molecu- 
lar clouds, or active Galactic nuclei (AGN). However, there are 
still a few TeV sources that could not be unambiguously associ- 
ated with any source de tected in lower energy bands (see, e.g., 
lAharonian et al. 2008a). In recent efforts to investigate these 
'dark' VHE emitters, multi-wavelength observations were con- 
ducted to probe the envir onment of these particl e accelerators 
(for recent work, see, e.g. . lAharonian et al E^bc). 

HESS J 1626-490 is another VHE gamma-ray source of un- 
known origin, which so far could not be identified with a 
source at lower wavelengths. This object, with an intrinsic ex- 
tension of ~5 arcmin (Gaussian FWHM), is located right on 
the Galactic plane (R.A.: 16 / *26 m 04*, Dec: -49°05'13") and 
was detected by H.E.S .S. with a peak significance of 7.5<x 
(lAharonian et al I l2QQ8ah . It is also stated that the source might 
be composed of two separate sources owing to the long tail 



extending towards the east. These authors measured a power- 
law spectrum with a photon index of 2.2 ± 0A stat ± 0.2 sys 
and a flux normalization of (4.9±0.9)xl0~ 12 cm~ 2 s" 1 TeV" 1 at 
ITeV at energies between 0.5 TeV and 40 TeV. A number 
of possible counterparts fro m other wavelengths are discussed 
by (lAharonia n et al ] 12008 A . such as the nearby faint exte nded 
ROSAT source 1RXS Jl 62504 -4909 18 dVoees et al.| [f999). the 
shell-type SNRG335.2+00.1 dWhiteoak & Greed Il996h . and 
despite their large offsets, the X-ray binaries (XRBs) 4U 1624- 
490 and IGR 16283-4838. Furthermore, the high-energy (HE) 
source 1FGL J1 626.0-4917c from the one-year Fermi source 
catalog (1FGL, lAbdo et aTl 120101) lies in close proximity to 
HESSJ1626-490. 

The two physical mechanisms thought to produce VHE 
gamma rays are inverse-Compton (IC) upscattering of low- 
energy photons by a population of relativistic electrons, as 
well as production and subsequent decay of 7r°s when high- 
energy hadrons interact with a dense medium. The latter, purely 
hadronic scenario becomes more viable in cases where the 
VHE emission correlates with the location of dense molecular 
clouds (e.g. towards the Galactic Center Ridge. lAharonian et al.l 
l2006bl) or when SNRs interact with the surrounding inter- 
stellar medium (ISM) as seen with H.E.S.S. from W28 and 
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other SNRs (lAharonian et al.ll2008d: iFeinstein etaDl2009l) and 
with VERITAS an d MAGIC from IC 443 (lAcciari et al.l 120091: 
I Albert et aDl20 07). In an IC scenario, X-ray emission coinciding 
with the VHE signal is expected due to synchrotron radiation 
from the same lepton population, whereas only a very low X-ray 
flux is predicted for a purely hadronic scenario arising from syn- 
chrotron cooling and/or nonthermal bremsstrahlung from sec- 
ond ary electrons prod uced in the decay of charged pions (see 
e.g. lGabicietaDl2Q09h . 

In this work we analyzed the data of an archival XMM- 
Newton observation to search for an X-ray counterpart of 
HESS J1626-490. Therefore, we classified the detected X-ray 
point sources in the vicinity of HESS J 1626-490 and searched 
for possible diffuse excess emissions above the expected 
Galactic background. Furthermore, we present 12 CO(J=1-0) 
molecular line survey data taken with the Nanten mm/sub-mm 
observatory to scan for molecular clouds and infrared data from 
the Spitzer GLIMPSE and MIPSGAL surveys to search for indi- 
cations of recent star-forming activity. 




i 16:26:00.0 
RA (J2000) 



2. X-ray data analysis and results 

We retrieved the EPIC (European Pho ton Imaging cam- 
era) data of an archival XMM-Newton dJansen et al.l 1200 ll) 
observation (ID: 0403280201), poi nted to the positio n of 
HESS J162 6-490. The EPIC- MOS (iTurner et al.l 120011) and 
EPIC-PN (IStruder et al.1 l200ll) data were analyzed with the 
XMM-Newton Science Analysis System (SAS) version 9.0.0, 
supported b y tools from th e FTOOLS package and XSPEC ver- 
sion 12.5.0 (lArnaudl 19961) for spectral modeling. For image pro- 
cessing, we used some tools from the CIAO 4.0 software pack- 
age. 

This observation was affected by long intervals of strong 
background flaring. To clean the data we applied a background 
good-time-interval (GTI) screening based on the full field-of- 
view (FoV) 7-15 keV light curve provided by the standard pro- 
cessing chain. We used thresholds of 8 cts/s for PN and 3 cts/s 
for MOS and the resulting net exposures are 4.9 ks for PN and 
13.2 ks for MOS, respectively. For all spectra and images pre- 
sented in the following sections, we selected good (FLAG==0) 
single and multiple events: PATTERN<4 (PN), PATTERNS 12 
(MOS). 



2.1. X-ray point sources 



Even though lAharonia n et al ] (I2008ah detected 
HESS J1626-490 with an intrinsic extension of -0.1°, we 
first characterized possible point-like X-ray counterparts. 
The SAS standard maximum likelihood technique for source 
detection was used in several energy bands: 0.2-0.5 keV, 
0.5-1.0 keV, 1.0-2.0 keV, 2.0-4.5 keV, 4.5-10.0 keV, and 0.5- 
10.0 keV. The resulting point source list indicates a detection 
threshold of 2x10" 14 erg cm -2 s" 1 . Within the 4<x contours 
of HESS J1626-490, we detected twelve X-ray point sources 
(see FigQ]). For five of these sources (1, 5, 7, 8, and 9) the 
photon statistics were sufficient to perform a spectral analysis. 
The remaining seven objects were classified according to their 
hardness ratio (HR), which we defined as the normalized differ- 
ence between high-energy (Hi: 2.0-10.0 keV) and low-energy 
(Low: 0.2-2.0 keV) counts: 



HR 



Hi - Low 
Hi + Low 



(1) 



Fig. 1. Combined color-coded image from the PN and MOS de- 
tectors. Energy intervals for the colors are 0.5-1.0keV (red), 1.0- 
2.0 keV (green), and 2.0-4.5 keV (blue). The red contours de- 
note the 7,6,5,4 o~ statistical significance levels of the H.E.S.S. 
detection of HESS J1626-490. Marked are point sources de- 
tected within the 4<x VHE contours (numbered 1 to 12 with in- 
creasing R.A.). The dashed ellipse (white) shows the lcr confi- 
dence region of the positional uncertainty of the Fermi source 
1FGLJ1626.0-4917C. 



The number of counts in the respective energy bands is the sum 
over all instruments where the source is detected. A correc- 
tion for differential vignetting was applied to the HRs. Similar 
to the Galactic plane X-ra y source populatio n studies con- 
ducted with XMM-N ewton (lHands et al.l 120041) and Chandra 
(lEbisawa et al.120 05). we classified sources as either soft (HR <- 
0.5), medium (-0.5< HR <0.5) or hard (HR >0.5). To identify 
near-infrared (NIR) counterparts, we correlated t he X-ray point- 
sourc e positions to the 2MASS all-sky catalog (iSkrutskie et al.l 
2006). Especially for regions on the Galactic plane where 
the interstellar absorption is very strong which is the case 
for HESS J1626-490 (N H ~2.22xl0 22 cnr 2 lDickev & Lockmanl 
Il990bl) . the NIR band is particularly suited to probing deep into 
the interstellar medium (ISM). We classified an X-ray source 
from this observation as NIR-identified when a 2MASS cata- 
log entry is found within 2 arcsec of the source position re- 
sulting from the SAS source detection algorithm. In the fol- 
lowing section we classify the detected point sources as either 
"soft Galactic" (e.g., X-ray active stars), "hard Galactic" (e.g., 
cataclysmic variables (CVs), X-ray binaries), or "hard Extra- 
galactic" (e.g., background AGN), based on the X-ray spectral 
characteristics and the NIR counterpart. These classifications 
are meant to reflect general characteristics such as hard/soft or 
absorbed/unabsorbed, which helps for deriving more accurate 
fluxes for sources where no spectral analysis was possible. 



2.1 .1 . Soft Galactic sources 

Owing to coronal activity and/or binary interaction, stars can 
exhibit strong thermal X-ray emission during a wide variety of 
evolutionary stages (for a com pilation of the 10 brightest X-ray 
stars detected by ROSAT, see lMakarovll2003l) . Therefore, such 
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Fig. 2. EPIC-PN and EPIC-MOS spectra of source No. 1 with a 

2 

3-temperature thermal plasma model (MEKAL, ^ = 1.07) fit. 

objects comprise the vast majority of soft Galactic sources. A 
typical physical model to describe their X-ray spectra is an ab- 
sorbed multi-temperature thermal plasma (MEKAL in XSPEC) 
with temperatures in the range of kT = 0.5-2.0 keV. Due to the 
lack of high interstellar absorption at these short distances and 
their intrinsically soft spectra, these objects would be classified 
as soft, according to the HR criterion, mentioned in the previ- 
ous section. Furthermore, we required an NIR counterpart for an 
identification as a "soft Galactic" source. According to these cri- 
teria, eight of the twelve point sources were classified as "soft 
Galactic", and a spectral analysis could be performed on three 
of them. We calculated the flux for those sources with no spec- 
trum from the 0.5-10.0 keV count-rate assuming an unabsorbed 
power-law spectrum with photon index 3.0. The X-ray and in- 
frared information for these objects are compiled in Table [T] 

Source No. 1 is by far the brightest X-ray point source 
(Fx,unabs = 8.0xl0" 12 erg cm" 2 s" 1 ) in the FoV and is spa- 
tially coincident with the active triplet system HD 147633 = 
1RXS J162517.7-490855(|M^Qv|2003|). This system contains 
a close binary of two G-type main sequence stars. The X-ray 
spectrum (Fig. [2]) can be well-fit with a 3-temperature thermal 
plasma emission model (see Table [T]). The three plasma temper- 
atures resulting from the fit are compa tible with tempera tures 
seen, e.g., from the Castor X-ray triplet (iGudel et al .1120011) . We 
therefore conclude that an identification of Source No. 1 with 
HD 147633 is very likely. 

2.1 .2. Hard Galactic and extra-galactic sources 

In contrast to the X-ray emission seen from nearby active stars, 
the spectra of background AGN and Galactic CVs are intrin- 
sically harder, and, in the case of AGN, heavily absorbed by 
the interstellar medium in the Galactic plane. According to 
the hardness-ratio classification scheme, these sources would 
be medium or hard. For two sources, not classified as "soft 
Galactic", an X-ray spectrum could be extracted (Nos. 7 
and 9). Both feature hard power-law spectra. Source No. 7 
shows a very high absorption column density (Nh = 8.8xl0 22 
cm -2 ), which is highe r than the Galactic level in that direction 
(N H = 2.1xl0 22 cm" 2 . iDickev & Lockmanlll990bh . so we clas- 
sify this object as "hard extra-galactic". In contrast, the spec- 
trum of the NIR identified source No. 9 is much less absorbed, 
indicating a location in the Galactic plane. We therefore clas- 
sify this source as "hard Galactic". The rest of the medium or 



hard sources were classified as "Galactic" if an NIR counterpart 
is present, otherwise as "extra-galactic". We calculated the flux 
for sources with no spectrum from the 0.5-10.0 keV count rate 
assuming an absorbed power-law spectrum with photon index 
2.0. For "extra-galactic" sources we assumed the total Galactic 
column density Nh = 2.1xl0 22 cm -2 , for "Galactic" sources the 
somewhat lower value Nh = 5xl0 21 cm" 2 . The results for all 
point sources are shown in Table [T] We note here that in the 
rare case of hard Galactic sources with no NIR counterpart, such 
as PWN, these objects would have been wrongly classified as 
"extra-galactic". 



2.2. Properties of the diffuse X-ray emission 

(lAharonian et al.1 l2008al) approximated the morphology of 
HESS J1626-490 by fitting a 2-D Gaussian to the VHE excess 
map. The resulting intrinsic (with the effects of the instrument 
point- spread-function removed) major and minor axes were 0.1 
deg and 0.07 deg, respectively, with a position angle of 3 deg 
west to north. An IC scenario would predict X-ray emission 
due to the synchrotron cooling of the same population of elec- 
trons. Likely candidates for such a scenario are PWN. For these 
sources the X-ray nebula is more compact than the VHE source 
since the high magnetic fields close to the pulsar lead to strong 
synchrotron cooling. In contrast, the VHE nebula might be dom- 
inated by a more extended population of electrons with lower 
energies and longer cooling timescales. However, as we did not 
see any obvious extended X-ray counterpart with XMM-Newton 
we characterized the diffuse X-ray emission seen from a region 
comparable to the VHE emission. 

Before the extraction of events from diffuse regions, we first 
removed all detected point sources from the dataset. For that pur- 
pose we mostly relied on the results from the SAS maximum- 
likelihood source detection algorithm. Exceptions were very 
bright sources (e.g. source No. 1), sources at large off-axis an- 
gles (e.g. source No. 5), and the out-of-time event column from 
source No. 1, where we modified or added the respective regions 
manually. All exclusion regions are shown in Fig. [3] 

To estimate of the non-X-ray background (NXB) we used 
a filter- wheel cl osed dataset, provide d by the EPIC background 
working group (ICarter & Readll2007h . To clean the background 
data of soft proton flares and of contributing observations with 
a general higher flux level, we applied the same GTI thresh- 
old as for the source data (see Sect. [2]). However, due to the 
non- simultaneity of the source and background observations, 
the particle-induced background level could still be different. 
At energies above ~12 keV the effective area of the X-ray tele- 
scope decreases significantly and the NXB becomes the domi- 
nant background component. Thus to improve the estimation of 
the NXB we derived a scaling factor based on the ratio of the 
count-rate in the 12-15 keV energ y bands between source an d 
background observations (similar to lDe Luca & Molendill2004l) . 

To check for possible hard diffuse X-ray excess emission, we 
extracted an image in the 3-7 keV energy band. We refilled ex- 
cluded regions using dmfilth and applied the adaptive smoothing 
algorithm c smooth (Fig. [3]). For the kernel size of the smoothing 
we required a minimum significance of 5<x. As background we 
subtracted the scaled NXB dataset as described in the previous 
paragraph. 
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Table 1. Classification of point sources 



Source 


R.A. 


Dec. 


N H 


HR^ 


kT/r (2) 


F (3) 
A,unaos 


J(4) 


2MASS 


class 


No. 






(xlO 21 cm- 2 ) 




(for kT: keV) 




(mag) 


name 




1 


16:25:17.5 


-49:08:53 


<0.08 


-0.91 


kT 2 = 0.20^ 
kT 2 = 0.63^ 
kT 3 = 2.0^g 


79.5 


6.117 


J16251768-4908524 


soft Galactic 


2 


16:25:19.6 


-49:13:33 


— 


-0.80 


— 


0.39 


13.458 


J16251980-4913341 


soft Galactic 


3 


16:25:45.0 


-49:04:30 


— 


-0.05 


— 


3.5 




— 


med. extra-galactic 


4 


16:25:55.4 


-49:12:56 


— 


-0.25 


— 


2.1 


12.895 


J16255550-4912576 


hard Galactic 


5 


16:26:08.0 


-49:00:09 


<1.6 


-0.72 


kTi = 

kT 2 = 3.5 + _\- 5 4 


3.1 


7.884 


J16260819-4900103 
- 


soft Galactic 


6 


16:26:09.2 


-49:05:21 




-0.82 




0.044 


12.42 


J16260907-4905195 


soft Galactic 


7 


16:26:17.2 


-49:09:03 




+0.92 


r = o.4i*g 


3.9 






hard extra-galactic 


8 


16:26:21.6 


-49:08:33 


<0.9 


-0.91 


kT = 0.84+^ 


1.9 


13.133 


J16262176-4908333 


soft Galactic 


9 


16:26:22.3 


-49:02:48 


6.5+5.0 


-0.16 


r = 2.0+0.7 


1.4 


11.304 


J16262255-4902485 


med. Galactic 


10 


16:26:38.1 


-48:59:29 




-0.98 




0.091 


13.47 


J16263792-4859277 


soft Galactic 


11 


16:26:47.2 


-49:00:22 




-0.39 




2.1 


12.726 


J16264736-4900230 


hard Galactic 


12 


16:26:57.4 


-49:10:30 




-0.73 




0.35 


11.704 


J16265766-4910312 


soft Galactic 



(1) Hardness ratio as defined in the text. (2) Temperature(s) or photon index resulting from a MEKAL or power-law model fit to the spectra; depending 
on the statistical quality of the thermal spectra, the number of temperature components varies. (3) Unabsorbed flux in the 0.5-10.0 keV band (xlO -13 
erg cm -2 s _1 ). The flux is derived from spectral fitting or, for fainter sources, by scaling an assumed spectrum with the count-rate depending on 
the source class (see text). (4) J magnitude from the 2MASS catalog in cases where a counterpart was found within 2 arcsec of the X-ray source. 




246.8 246.6 246.4 246.2 

Right ascension 

Fig. 3. Adaptively smoothed EPIC-PN image of diffuse X-ray 
flux in the 3-7 keV energy band with a linear color scale (ar- 
bitrary units). The instrument FoV is shown as a dashed circle 
(green). The large ellipses (black) denote the extraction regions 
for HESS J 1626-490 and the background test, respectively. The 
large circle (magenta) to the east gives the approximate exten- 
sion of the X-ray scattering halo of 4U 1624-490. Excluded re- 
gions are shown as crossed-out areas (yellow regions crossed 
with red lines). 







fefet, 

T+Tf 




+■ 


WW 





















1 2 5 

Channel Energy (keV) 



Fig. 4. NXB subtracted spectra from the elliptical region cen- 
tered at HESS J1626-490 (black) and a region in the south of 
the FoV outside the 4 cr contours of HESS J1626-490 (red). The 
stepped lines show the best-fit of a 2-temperature NEI model. 



2.2.1 . Upper limit for extended X-ray emission from 
HESS J1 626-490 

The diffuse image (Fig. [3]) shows no significant excess emission 
above the flat Galactic diffuse component apart from the area at 
the eastern edge of the FoV, which is discussed in the next sec- 
tion. To derive an upper limit for the total X-ray flux possibly 
connected to HESS J1626-490 (yellow ellipse in Fig. [3]) we ex- 
tracted a spectrum from that area and compare it to a spectrum 
extracted from the southern part of the FoV outside the 4 cr con- 
tours of the VHE detection. For the spectral analysis we used the 
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same exclusion regions as for the extraction of the diffuse image 
(see Sect. CHI). 

To accurately account for the energy-dependent variation of 
the effective area over th ese large regions we applied the weight- 
ins method described bv lArnaud et"atl(l200lh to extract the spec- 
tra. In accordance with this method we used the on-axis energy 
response (RMF) and ancillary response files (ARF). As back- 
ground we subtracted spectra from the identical regions in the 
NXB dataset using the weighting method as well. Both NXB 
subtracted spectra are shown in Fig.|4l 

Since HESS J1626-490 is right on the Galactic plane, 
we assumed that the observed emission from both regions 
is of diffuse Galactic origin. Therefore, we adopted the two- 
temperature non-equ ilibrium ionization model (2-T NEI) from 
lEbisawa et al] (120051 hencefor th E05). These authors obser ved a 
typical Galactic plane region (iDickey & Lockmanlll990bl) with 
Chandra and analyzed the Galactic diffuse emission in great de- 
tail. We fixed most of the model parameters at the values from 
E05 except for the norm and column densities of both temper- 
ature components and the Si abundance of the hard component, 
which would have been significantly overestimated otherwise. 
We fit both spectra in parallel with this model (Fig. [4]), and we 
only allowed the normalizations to vary separately. The intrinsic 
surface fluxes were derived by dividing the unabsorbed fluxes by 
the effective extraction region, which is the region on the detec- 
tor minus excluded areas, bad pixels/columns, and CCD gaps. 
For both regions the details of the fitting parameters are listed in 
Tabled 

The total surface fluxes for both regions agree within 1.5 <x, 
with the flux measured from the test region slightly larger. 
Therefore, we did not detect any significant X-ray excess emis- 
sion from the direction of HESS Jl 626-490 with respect to the 
test region. Assuming that a possible X-ray signal associated to 
HESS J 1626-490 would originate in a region with the same ex- 
tension as the VHE signal and that an excess flux greater than 
4<x above the Galactic diffuse emission would be detectable, 
we derived an upper limit for the X-ray flux in the 1-10 keV 
band of F X5ex cess < 4.85xl0" 12 erg cm -2 s" 1 . Given the observed 
background count rate from the ellipse and assuming an ab- 
sorbed power-law spectral shape with index -2 and column den- 
sity lxlO 22 cm -2 , the formal 4<x confidence level flux upper 
limit (1-10 keV) would be ~lxl0" 12 erg cm -2 s" 1 , somewhat 
lower than the above value. However, we used the first value as 
we deem it more realistic for the actual source analysis tech- 
niques. 

Figure \5\ shows the VHE spectral energy distribution (SED) 
of HESS J1626-490, together with this upper limit, assuming 
a power-law spectrum with index -2 for the X-ray spectrum. 
In addition, the spectral uncert ainty band of the Fermi source 
lFGLJ1626.0-4917cis shown (lAbdo et al.ll20Toh . 

2.2.2. Contamination from the dust-scattering halo of the 
low-mass X-ray binary 4U 1624-490 

As seen in the smoothed image of diffuse X-ray emission [3l the 
flux level increases towards the eastern edge of the FoV. The dip- 
ping low-mass X-ray binary 4U 1624-490 is located -270 arcsec 
outside the XMM-Newton FoV, exactly to the east. This source 
features an exten ded dust- scattering halo, which was detected 
with BeppoSAX fe alucinsk a-Church et al. 1 12000b and Chandra 
dXiang et al ] l2007l) 7 According to these authors the halo shows 
a hard, highly absorbed spectrum (Nr = 8-9xl0 22 cm -2 ) and 
should reach well inside the XMM-Newton FoV of this observa- 
tion as denoted by the eastern yellow circular region in Fig. [3] 
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Fig. 5. Spectral energy distribution of HESS J1626-490, show- 
ing the H.E.S.S. spect rum as data points with errors (red) 
(lAharonian et al.ll2008al) togeth er with the spectra l uncertainty 
band of 1FGL J1626.0-4917c dAbdo et alJl2010b . The XMM- 
Newton upper limit from this work is denoted by a horizontal 
line with an arrow (black). 



An absorbed power-law fit to the spectrum extracted from the 
eastern region in the XMM-Newton observation gives a column 
density of Nh = (8.0±2.5)xl0 22 cm" 2 and a photon index of 
T = 2.3+0.7. These results, together with the apparent circular 
shape of the excess, confirm the identification of this emission 
with the dust- scattering halo of 4U 1624-490. The effects of this 
contamination are very limited to the eastern edge of the FoV 
and should therefore not influence much the measurements of 
the diffuse emission from the direction of HESS Jl 626-490 and 
the test region. 



3. Nanten 12 CO(J=1-0) data 

To search for molecular clouds, spatially and morphologically 
coincident with HESS J1626-490, we analyzed 12 CO(J=1-0) 
molecular line observations performed by the 4 m mm/sub-mm 
Nanten telescope, locate d at Las Campanas Observatory, Chile 
( Mi zuno & Fukuill2004l) . These data were tak en as part of the 
Nante n Galactic Plane Survey (1999 to 2003. iMatsunaga et all 
l200ll and references therein). The survey grid spacing was 4 ar- 
cmin for the region around HESS J 1626-490. For this work the 
local standard of rest velocity (vlsr) range -240 to +100 km s" 1 
was searched. 

Figure [6] (Left) shows the Nanten 12 CO(J=1-0) image in- 
tegrated over the vlsr range -31 to -18 kms" 1 . In this inter- 
val we found a 12 CO feature partially overlapping with the 
VHE emission . Acco rding to the Galactic rotational model of 
Bra nd & Blitzl {[993), this vlsr range corresponds to a kine- 
matic distance of 2.2 to 1.5 kpc. As can be seen in the Nanten 
position- velocity plot (Fig.\6\(Right)), this CO cloud seems to be 
connected to a much larger system of clouds at velocities -32 to 
-50 kms" 1 . The 12 CO emission at these more distant velocities 
is much more extended, and we did not find any apparent feature 
particularly matching the VHE morphology or coinciding with 
the SNRG335.2+00.1. 

Using the relation between the hydrogen column 
density (N H ) and the 12 CO(J=1-0 ) intensity W( 12 CO) , 
N H = 1.5 x 10 20 [W( 12 CO)/Kkm/s] dStrong et al.l 12004 . 
we estimate the total mass of this cloud at 1.8xl0 4 M o for d 
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Table 2. Fit results for a 2-T NEI model 



Parameter 



HESS J 1 626-490 Southern test region 



Soft Component 



kT (keV) 
log(n e t) (cm~ 3 s) 
Abundance (except Ne, Mg, Si) 
Ne abundance 
Mg abundance 
Si abundance 
N H (xlO 22 cm" 2 ) 
Intrinsic surface flux ( * } (xlO 



0.59 (frozen) 
11.8 (frozen) 
0.044 (frozen) 
0.30 (frozen) 
0.14 (frozen) 
0.25 (frozen) 
0.23+0.09 



7 erg cm 2 s 



sr 1 ) 



1.4+0.4 



1.8+0.5 



Hard Component 



kT (keV) 
log(n e t) (cm~ 3 s) 
Abundance (except Fe) 
Fe abundance 
N H (xlO 22 cm" 2 ) 

Intrinsic surface flux ( * ) (xlO -6 erg cm -2 s~ 

x 2 /y 



sr 1 ) 



2.7+0.2 



5.0 (frozen) 
10.6 (frozen) 
0.17 (frozen) 
0.9 (frozen) 

o 17+O.31 
D ' 1 -0.27 

102/74 



3.6+0.3 



w In the 0.7-10.0 keV energy band 




Fig. 6. Left: Nanten 12 CO(J=1-0) image of the region around HESS J1626-490 (linear scale in Kkms -1 ) integrated over the vlsr 
range -31 to -18 km s" 1 . Overlaid are the contours of the VHE emission (blue) and of the H I cloud discussed in Sect. |4] (green). The 
dashed circle (light blue) denotes the position and extension of SNR G335.2+00.1. Right: Nanten 12 CO(J=1-0) Galactic longitude- 
velocity plot (linear scale in Kdeg) integrated over the Galactic latitude range -0.33 to 0.14 deg. The 12 CO cloud shown in the left 
image is marked by a dashed circle (green). 



= 1.8 kpc within an elliptical region centered at / = 334.78, 
b = 0.00 with dimensions 0.26x0.30 deg. The corresponding 
average density is 2.1xl0 2 cm -3 . 

4. HI and infrared data 

Figure [8] shows HI images of the region around 
HESSJ1626-490 from the Southern Galactic Plane Survey 
(SGPS iMcClu re-Griffit hs et al.ll2005h integrated over the v L sr 
ranges -23 to -18 kms" 1 {Left) and -31 to -23 kms" 1 (Right), 
respectively. A local H I density depression is seen in the center 



of the left image at an angular separation of -21' from the CO 
molecular cloud and HESS J1626-490. It is striking that this 
feature is consistent in positi on as well as i n angular extensio n 
with the SNRG335.2+00.1 (iGreenl 120091: iReach et alJ l2006h . 
Furthermore, both images show a region of increased gas 
density in spatial coincidence with the VHE contours, which is 
more pronounced in the -31 to -23 kms -1 velocity range. Two 
dips, most likely due to foreground absorption, are marked by 
arrows (labeled Dl and D2). 

We estimated the mass of the dense HI region coincid- 
ing with the VHE and 12 CO features using the signal within 
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SNR G335.2+00.1 



E „ 



o 

<D 
> 



335.6 



335.4 335.2 335 334.8 
Galactic Longitude (deg) 




334.6 



Fig. 7. H I SGPS Galactic longitude-velocity plot (linear scale 
in Kdeg) integrated over the Galactic latitude range -0.11 to 
0.24 deg. The position and extension of the 12 CO cloud is de- 
noted by a dashed circle (black). The double arrow (green) 
shows the Galactic latitude and extension of SNR G335.2+00.1. 

the same elliptical region as for the CO cloud (Sect. [3]). Using 
the relation betwee n H I in tensity and column density from 
iDickev & Lockmanl (Il990al) (X=1.8xl0 18 cm" 2 K" 1 km/s) we 
estimated the mass of the cloud as 4.9xl0 3 M© with an aver- 
age density of 60.1 cm" 3 . However, we note here that due to the 
two absorption dips in the region this value should be seen only 
as a lower limit of the actual mass of the H I cloud. 

Figure|7]shows the vlsr profile for this region integrated over 
the Galactic latitude range -0.1 1 to 0.24 deg, which is the exten- 
sion of both the SNRG335.2+00.1 and the VHE emission of 
HESS J1626-490. Both, the region of low density and the HI 
cloud are clearly visible. The latter is in good agreement with 
the position of the CO cloud (Sect. [3]), as indicated by the dashed 
ellipse. 

Figure [9] shows a three-color (rgb= 24/8/3. 6 fim) image from 
the Spitzer GLIMPSE and MIPSGAL surveys in log scale 
(MJy/sr units) with lower limit clipping at (rgb= 3 0/60/5 MJy/sr). 
Apart from the large-scale diffuse emission in the red and green 
bands along the Galactic plane, there are two prominent HII 
regions within the 4<x H.E.S.S. contours. These two HII re- 
gions are coincident with the absorption dips seen in H I (Fig. [8]), 
which is a com mon feature of many H I I regions in the Galactic 
plane (s ee, e.g.jKuchar & Banialll990l) . Rl is listed in the cat- 
alog of iRusse il (2003) (No. 427) with a systemic distance of 
2.4+0.3 kpc. Even though R2 is not listed, judging from the sim- 
ilarity of the two absorption dips, the distances to both H II re- 
gions should be quite similar. Both HII regions are also clearly 
seen in th e Molon^lo (|Green et al 1 fl999h 843 MHz radio image 
shown by (Ah aronian et al . 2008a, Fig. 2). The differing distance 
estimates for the H II regions and the H I/CO-complex might in- 
dicate systematic differences between both methods. In particu- 
lar this could mean that a large fraction of the mass of the H I 
cloud is located at a distance more than 1.8 kpc. 

5. Discussion 

In this work we searched for X-ray, sub-millimeter, and in- 
frared counterparts of the unidentified VHE gamma ray source 
HESS Jl 626-490 to explore the nature of the particle acceler- 
ation process responsible for the observed emission. With a lu- 



minosity of 2xl0 33 erg s _1 x(d/kpc) 2 , HESS J1626-490 is quite 
typical of Galactic TeV sources detected by H.E.S.S., and we 
discuss both leptonic and hadronic emission scenarios in this 
section. 

According to lReimer et al.l (120061) . close binary systems con- 
sisting of Wolf-Rayet (WR) and/or OB super-giant stars are ex- 
pected to effectively produce VHE gamma-rays in the colliding 
wind zone where strong shocks are present and where particle 
acceleration to multi TeV energies could take place. The nearby 
triplet system HD 147633 (source No. 1, Sect. 12. Lit is a strong 
X-ray emitter. Its spectrum can be well fit by a thermal model 
with temperatures typical of X-ray binaries showing coronal ac- 
tivity. This system is composed of a close binary of two G-type 
main sequence stars a nd a tertiary star orbiting at a greater dis- 
tance (Maka rovl2003l) and can therefore most likely be ruled out 
as an efficient source of highly relativistic particles, which could 
give rise to VHE gamma-ray emission. Apart from that, the ex- 
tended morphology of HESS J 1626-490 makes an identification 
with an X-ray point- source detected by XMM -Newton rather un- 
likely. 

In the case of a leptonic scenario where low-energy pho- 
tons are up- scattered by relativistic electrons via the IC pro- 
cess, X-ray emission is expected to accompany the VHE sig- 
nal arising from synchr otron cooling of the s ame population 
of high-energy electrons (I Aharonian et al.|[l997l) . Two examples 
for such sources d etected by H.E.S.S. are t he pulsar wind neb- 
ulae MSH 15-52 (I Aharonian et al.l l2005bl: ISchock et al.l l2010h 
and HESS J1825-13 (lAharonian et al.l l2005cl). Even from a first 
glance at the flux levels in X-rays (~xl0 -13 erg cm" 2 s" 1 for 
point sources and ~xl0" 12 erg cm" 2 s" 1 for diffuse emission) 
and at VHE energies (~xl0" n erg cm" 2 s" 1 ), a purely leptonic 
emissi on scenario seems unlikely. Following lAharonian et al.l 
(fl997h and assuming a typical Galactic magnetic field strength 
of 3xl0" 6 G, the corresponding synchrotron photon energy for 
E=l TeV IC photons is 6=0.02 keV. Assuming the same photon 
index for the synchrotron counterpart as measured from the VHE 
source (T = 2.18) and extrapolating the flux to the 0.5 to 10 keV 
band, we estimate an integrated source flux of -1.1x10" 11 erg 
cm" 2 s" 1 . This flux is a factor of ~25 more than what we measure 
from any point source, apart from X-ray source No. 1 , which was 
discussed in the previous paragraph. This makes an identification 
of HESS J 1626-490 with any of these sources unlikely. Our up- 
per limit for diffuse X-ray emission (Sect. lZXTt is a factor of ~2 
lower than the expected value for HESS Jl 626-490 in a purely 
leptonic model. 

We observed infrared emission from the field around 
HESS Jl 626-490 (see Sect|4]) which could provide an additional 
target radiation field for the IC process. This would lower the 
expected synchrotron flux in the X-ray band. On the other hand, 
these estimates are based on the lowest possible magnetic field 
found in the Galactic ISM. Magnetic fields in the vicinity of en- 
ergetic pulsars can be a factor of 10-100 larger, which would 
increase the expected synchrotron flux accordingly. We there- 
fore conclude that a purely leptonic emission scenario is rather 
unlikely for HESS J1626-490. 

Not detecting any X-ray source fulfilling the energetic re- 
quirements for a purely leptonic scenario favors a hadronic emis- 
sion process such as dense clouds in the vicinity of a cosmic par- 
ticle accelerator. Such a scenario will be discussed in the remain- 
ing part of this section. As already mentioned in Sect. [T] dense 
molecular clouds are established VHE gamma-ray emitters be- 
cause they provide target material in regions of high cosmic- 
ray densities. Using data from the Nanten 12 CO(J=1-0) Galactic 
plane survey, we detected a molecular cloud that is morpholog- 
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335.5 335 334.5 335.5 335 334.5 

Galactic longitude (deg) Galactic longitude (deg) 



Fig. 8. HI SGPS images of the region around HESS J1626-490 (linear scale in Kkms" 1 ) integrated over the vlsr ranges -23 to 
-18 kms" 1 (Left) and -31 to -23 kms" 1 (Right), respectively. The dashed circle (white) marks the position and radio extension of 
SNRG335.2+00.1. Also shown are the H.E.S.S. VHE contours of HESS J1626-490 (green). The arrows (white, labeled Dl and 
D2) mark two dips, probably due to foreground absorption. 





Fig. 9. Three-color (rgb=24/8/3.6 /zm) image from the Spitzer 
GLIMPSE and MIPSGAL surveys in log scale (MJy/sr units) 
with lower limit clipping at (rgb= 3 0/60/5 MJy/sr) in Galactic 
coordinates. Shown are the VHE contours of HESS J 1626-490 
(white) and two HII regions (red, labeled Rl and R2). 



ically consistent with HESS J1626-490. This object is located 
at a ki nematic distance of -1.8 kpc. Following I Aharonian et al.l 
(11994 Eqs. 2 and 3), we estimated the required gas density 
to produce the observed VHE y-ray signal (F r (>0.6TeV) = 
7.5x10" 12 ph cm -2 s" 1 and Y = 2.2) as n « 10 cm -3 assuming 
a cosmic ray production efficiency of 6 = 0.1 and a distance of 
d = 1.8 kpc. This value is an order of magnitude lower than the 
measured 12 CO mean density. Thus, this environment would be 
easily suited to providing the observed VHE y-ray flux. 



Following Ah aroma 3 (119911) (Eq. 10), assuming that the 
observed molecular cloud only provides a 'passive' target for 
cosmic-rays originating in a nearby acceleration site, the cosmic- 
ray enhancement factor for this source would be kcR = 333 
for 10 TeV protons. As seen in Fig. [8] (Right) and Fig. [7] 
there is a region of high-density HI gas spatially coincident 
with HESSJ1626-490 located at v LSR ~ -25 kms" 1 . Taking 
the additional mass of this H I cloud into account as additional 
target material kcR would be down to -255 which is an or- 
der of magnitude l arger than seen from o ther SNR-ISM in- 
ter actions (see, e.g. jAharonian et al.ll2008d) . However, follow- 
ing [Ah^onian^Atoxao dl99g) (Fig. lb), even cosmic ray en- 
hancements of kcR > 10 3 are to be expected from molecu- 
lar clouds at distances <10 pc (as in our case, see next para- 
graph) from the SNR shell and for diffusion coefficients of 
£>(10GeV)~10 26 cm 2 s" 1 if the age of the SNR does not exceed 
10 4 to 10 5 years. 

In particular in the left image of Fig.[8j there appears to be a 
region of higher concentration of H I gas towards the eastern tail 
of the VHE emission, which in turn could arise from a second 
VHE source (as already mentioned in Sect.[T]) . The 12 CO emis- 
sion that traces the denser regions of the cloud does not show 
two distinct features in a similar way. As this issue cannot be 
clearly solved based on the current available VHE data, we did 
not investigate any subregions and only considered the area as a 
whole. 

Now assuming that this 12 CO /HI cloud is indeed the source 
of the observed VHE gamma-rays, a nearby accelerator would 
be needed. The presence of two prominent H II regions at dis- 
tances consistent with the CO /HI clouds (see Sect.© indicates 
massive star formation, hence potential PWNe and SNRs in this 
area. The nearby density depression seen in H I (see Sect. |4]) 
might indicate the presence of a recent catastrophic event, such 
as an SNR, giving rise to strong shocks that would have blown 
the neutral gas out. At d = 1.8 kpc the edge of this region would 
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be at a distance of 8.1 pc from the 12 CO/HI cloud. It is strik- 
ing that this H I feature is consistent in bot h position an d angu- 
lar extension w ith the SNRG3 35. 2+0 0.1 (lGreenll2009l) . Using 
the X-D relation Guseinov et al ] <l2004h estimated the distance to 
SNRG335.2+00.1 as d = 5 kpc. which conflicts with the iden- 
tification of SNRG335.2+00.1 with the HI depression seen at 
a distance of 1.8 kpc. However, the scatter in E-D distances for 
individual sources is quite large (lGreenl l2005) making a positive 
identification of SNRG335. 2+00.1 and the HI depression seen 
at a kinematic distance of 1.8 kpc very well possible. 

We note here that, apart from possible systematic uncertain- 
ties in kinematic distance estimates, the projected distance of 
8.1 pc should be seen only as a lower limit because the phys- 
ical separation of accelerator and target could be larger if both 
objects are not at the same distance to the observer. A hint that 
this could be the case for this source is that the center of grav- 
ity of the CO cloud is slightly shifted in velocity with respect to 
the H I density depression. A larger distance between the accel- 
erator, most likely SNRG335. 2+00.1, and the molecular cloud 
would lower the expected cosmic ray enhancement. 

Even though the angular separation of the Fermi source 
1FGL J1626.0-4917c and HESS J1626-490 is relatively small, 
these sources are not consistent within the 3<x significance level 
of their spatial uncertainties (see Fig. Q]). Furthermore, we did 
not find any feature in X-rays, CO, or HI that could explain 
a displacement between HE and VHE gamma-ray emission. 
1FGL J1626.0-4917c shows signs of a slight variab ility. Based 
on th e variability index listed in the 1FGL catalog (lAbdo et al] 
120101) . the chance of the source being a constant emitter can be 
estimated to < 1 1 %, which might indicate an extragalactic origin. 
Thus it remains unclear at this point whether these two sources 
are physically connected. 

6. Outlook 

With the improved angular resolution of next-generation 
Cherenkov Telescope arrays, it will be possible to study the mor- 
phology of the VHE signal and its connection to molecular gas 
in much more detail. This is crucial for constraining the emission 
region of VHE y-rays, which strongly influences any estima- 
tion of target mass and density. For the determination of diffuse 
X-ray fluxes and upper limits, a large effective area combined 
with a good angular resolution to effectively exclude point-like 
sources is needed. Therefore, a much deeper exposure (-100 ks) 
with XMM-Newton would make it possible to significantly im- 
prove the constraints on lep tonic scenarios. However, for a more 
likely hadronic scenario, the flux expectations in the keV band 
are very low, which could make approval of such a proposal very 
unlikely. To study the particle propagation in more detail, dedi- 
cated radio polarization observations in the MHz to GHz regime 
would help determine the direction and level of turbulence of 
the magnetic field in the SNR. VHE gamma-ray emission from 
molecular clouds may also exhibit localized peaks at their cor e 
regions with densities >10 4 cm" 3 dGabici & Aharonianll2007l) . 
To trace these regions molecular lines of NH3, CS, or SiO can 
be used. Respective follow-up observations at mm and sub-mm 
wavelengths could be crucial for characterizing the environment 
of particle interaction. 

Acknowledgements. The XMM-Newton project is supported by the 
Bundesministerium fur Wirtschaft und Technologie/Deutsches Zentrum 
fur Luft- und Raumfahrt (BMWI/DLR, FKZ 50 OX 0001) and the Max-Planck 
Society. This publication makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California Institute of Technology, 



funded by the National Aeronautics and Space Administration and the National 
Science Foundation. The NANTEN project is financially supported from 
JSPS (Japan Society for the Promotion of Science) Core-to-Core Program, 
MEXT Grant-in-Aid for Scientific Research on Priority Areas, and SORST-JST 
(Solution Oriented Research for Science and Technology: Japan Science and 
Technology Agency). 



References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJS, 188, 405 
Acciari, V. A., Aliu, E., Arlen, T., et al. 2009, ApJ, 698, L133 
Aharonian, E, Akhperjanian, A. G., Aye, K.-M., et al. 2005a, Science, 307, 1938 
Aharonian, E, Akhperjanian, A. G., Aye, K.-M., et al. 2005b, A&A, 435, L17 
Aharonian, E, Akhperjanian, A. G., Barres de Almeida, U., et al. 2008a, A&A, 
477, 353 

Aharonian, E, Akhperjanian, A. G., Barres de Almeida, U., et al. 2008b, A&A, 
483, 509 

Aharonian, E, Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2008c, A&A, 481, 
401 

Aharonian, E, Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006a, ApJ, 636, 
777 

Aharonian, E, Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006b, Nature, 
439, 695 

Aharonian, F. A. 1991, Ap&SS, 180, 305 

Aharonian, F. A., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2005c, A&A, 
442, L25 

Aharonian, F. A. & Atoyan, A. M. 1996, A&A, 309, 917 

Aharonian, F. A., Atoyan, A. M., & Kifune, T. 1997, MNRAS, 291, 162 

Aharonian, F. A., Drury, L. O., & Voelk, H. J. 1994, A&A, 285, 645 

Albert, J., Aliu, E., Anderhub, H., et al. 2007, ApJ, 664, L87 

Arnaud, K. A. 1996, in ASP Conf. Ser. 101: Astronomical Data Analysis 

Software and Systems V, 17 
Arnaud, M., Neumann, D. M., Aghanim, N., et al. 2001, A&A, 365, L80 
Baluciriska-Church, M., Humphrey, P. J., Church, M. J., & Parmar, A. N. 2000, 

A&A, 360, 583 
Brand, J. & Blitz, L. 1993, A&A, 275, 67 
Carter, J. A. & Read, A. M. 2007, A&A, 464, 1155 
De Luca, A. & Molendi, S. 2004, A&A, 419, 837 
Dickey, J. M. & Lockman, F. J. 1990a, ARA&A, 28, 215 
Dickey, J. M. & Lockman, F. J. 1990b, Ann. Rev. Astron. Astrophys., 28, 215 
Ebisawa, K., Tsujimoto, M., Paizis, A., et al. 2005, ApJ, 635, 214 
Feinstein, E, Fiasson, A., Gallant, Y., et al. 2009, in American Institute 

of Physics Conference Series, Vol. 1112, American Institute of Physics 

Conference Series, ed. D. Bastieri & R. Rando, 54-62 
Gabici, S. & Aharonian, F. A. 2007, Ap&SS, 309, 465 
Gabici, S., Aharonian, F. A., & Casanova, S. 2009, MNRAS, 396, 1629 
Green, A. J., Cram, L. E., Large, M. I., & Ye, T. 1999, ApJS, 122, 207 
Green, D. A. 2005, Mem. Soc. Astron. Italiana, 76, 534 
Green, D. A. 2009, Bulletin of the Astronomical Society of India, 37, 45 
Gudel, M., Audard, M., Magee, H., et al. 2001, A&A, 365, L344 
Guseinov, O. H., Ankay, A., & Tagieva, S. O. 2004, Serbian Astronomical 

Journal, 169, 65 

Hands, A. D. P., Warwick, R. S., Watson, M. G., & Helfand, D. J. 2004, MNRAS, 
351,31 

Jansen, E, Lumb, D., Altieri, B., et al. 2001, A&A, 365, LI 
Kuchar, T. A. & Bania, T. M. 1990, ApJ, 352, 192 
Makarov, V. V. 2003, AJ, 126, 1996 

Matsunaga, K., Mizuno, N., Moriguchi, Y, et al. 2001, PASJ, 53, 1003 
McClure-Griffiths, N. M., Dickey, J. M., Gaensler, B. M., et al. 2005, ApJS, 158, 
178 

Mizuno, A. & Fukui, Y. 2004, in Astronomical Society of the Pacific Conference 
Series, Vol. 317, Milky Way Surveys: The Structure and Evolution of our 
Galaxy, ed. D. Clemens, R. Shah, & T. Brainerd, 59-+ 

Reach, W. T., Rho, J., Tappe, A., et al. 2006, AJ, 131, 1479 

Reimer, A., Pohl, M., & Reimer, O. 2006, ApJ, 644, 1118 

Russeil, D. 2003, A&A, 397, 133 

Schock, F. M., Biisching, I., de Jager, O. C, Eger, P., & Vorster, M. J. 2010, 

A&A, 515, A109+ 
Skrutskie, M. E, Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 
Strong, A. W., Moskalenko, I. V, Reimer, O., Digel, S., & Diehl, R. 2004, A&A, 

422, L47 

Striider, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 
Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 
Voges, W., Aschenbach, B., Boiler, T., et al. 1999, A&A, 349, 389 
Whiteoak, J. B. Z. & Green, A. J. 1996, A&AS, 118, 329 
Xiang, J., Lee, J. C, & Nowak, M. A. 2007, ApJ, 660, 1309 



